Band structure approach to the resonant x-ray scattering 
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We study the resonance behaviour of the forbidden 600 and 222 x-ray Bragg peaks in Ge using 
LDA band structure methods. These Bragg peaks remain forbidden in the resonant dipole scat- 
tering approximation even taking into account the non local nature of the band states. However 
they become allowed at resonance if the eigenstates of the unoccupied conduction band involve a 
hybridization of p like and d like atomic states. We show that the energy dependence of the resonant 
behaviour, including the phase of the scattering, is a direct measure of this p-d hybridization. and 
obtain quantitative agreement with experiment. A simple physical picture involving a product of 
dipole and quadrupolar transition matrix elements explains this behaviour and shows that it should 
be generally true for cases where the resonating atom is not at an inversion center. This has strong 
implications for the description of the resonance behavior of x-ray scattering in materials where the 
resonant atom is not at an inversion center such as V2O3 and in ferro and antiferro electric and 
piezo electric materials in general. 
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In a recent paper T.L. Lee et al g] describe a detailed 
study of the energy dependence of the 600 reflection in 
the vicinity of the Ge K-edge. They demonstrate a strong 
600 reflection close to the Ge K-edge with an energy de- 
pendence quite different from the absorption edge struc- 
ture in that the 600 intensity drops of rapidly at energies 
just beyond the edge. The intensity of this reflection at 
the edge corresponds to a structure factor of 0.45 elec- 
trons. The fact that one sees this reflection at all as 
well as the large intensity and the peculiar energy depen- 
dence is difficult to explain. This because, even if we take 
into account the non spherical nature of the charge dis- 
tribution and put electron charge in the bonds between 
the atoms from which the x-rays, scatter the 600 reflec- 
tion remains forbidden. Also the usual dipole-dipole or 
quadrupole-quadrupole like resonances as discussed for 
the manganites |, |, |, H and V2O3 [|, |, i cannot 
contribute to this reflection as discussed below. Recently 
Templeton and Templeton ^J suggested that the lowest 
order process which will allow such a reflection is one 
involving a product of a dipole and a quadrupole like 
matrix element for the scattering. In this paper we show 
that band structure calculations provide a natural and 
detailed explanation for the observed behaviour and that 
the resonance behaviour is a direct measure of the energy 
dependence of the p-d hybridization in the conduction 
band of Ge. This is, in a way equivalent to Templeton's 
suggestion but here is caused by a natural p-d mixing due 
to a lack of inversion symmetry at the Ge atomic sites. 

We also show that the 222 reflection becomes allowed 
in this fashion but this is not as spectacular since the 
222 reflection is also weakly allowed by simple Thomson 
scattering of the charge situated in the bonds between the 



atoms [|o| |ll|, |2|. Non the less the energy dependence 
of the phase of this extra scattering channel causes a de- 
crease in the net scattering close to the edge as observed 
and we predict an increase at energies about 30 eV above 
the edge because of a change in the phase of the p-d mix- 
ing at these energies. Although the case of Ge may, for 
some, not be that interesting it serves as a demonstra- 
tion of the effect and the success of theory to explain it. 
We suggest that this effect will play an important role in 
the resonant behaviour of recently discussed controversial 
materials like V2O3, Fe2 03 and many more vanadates , 
titanates and manganites which often also lack inversion 
symmetry at the resonating ion. In fact quite generally 
this effect could be very important in ferro and antiferro 
electric as well as piezoelectric materials. This may also 
provide a simple way to actually distinguish between ma- 
terials with and without an inversion center which is an 
important problem in crystallography. 

Ge has the diamond structure which contains 8 atoms 
per unit cell as depicted in Fig|| showing also the loca- 
tion of these atoms in a simple cubic notation. We notice 
immediately that if all the Ge atoms are equivalent then 
reflections like the 600 and 222 reflections will be for- 
bidden. Even if we take into account that the charge 
density around each atom is not spherical but somewhat 
enhanced along the bond directions we still will find zero 
intensity for the 600 reflection because of the space group 
symmetry. However in this case the 222 reflection can 
become weakly allowed in a purely Thomson scattering 
mechanism. Looking at Fig.n^more carefully however and 
at the Ge atoms connected by a bond it is easy to see that 
in fact Gei and Ge2 are not really equivalent in that the 
charge distribution around atom 1 is the, by 180 degrees, 




FIG. 1: The diamond crystal structure. 



rotated image of that around atom 2. In Thomson scat- 
tering this difference will not be noticed however. This 
difference between Gei and Ge2 is caused by a change 
in sign of the hybridization between s and p and the p 
and d wave functions forming the eigenstates of the oc- 
cupied valence and unoccupied conduction band states. 
For example if we take the Gei-Ge2 direction to be the z 
direction then a wave function like as + Ppz + 7^322-^2 
on Gei would have charge pile up in the Gei-Ge2 bond 
but for Ge2 this combination at the same energy would 
be as — j3pz + 'yd^z'^^r'^. Notice the sign change of the 
p contribution because of the ungerade character of its 
wave function. We note that s, p, d can mix in this struc- 
ture because of the lack of inversion symmetry about the 
Ge atoms. 

We will not see this difference in an experiment which 
is sensitive only to the charge density as in Thomson 
scattering. However in a resonant scattering process the 
anomalous form factor is given by 
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Where e' and e are the polarization of the scattered and 
incident radiation respectively and K' and K the respec- 
tive wave vectors. For the K-edge studies a refers to a Is 
orbital and 5 to a band eigenstate. 

Following Templcton we expand the exponential as 
1 + iK ■ f resulting in 3 kinds of contributions: dipole- 
dipole, quadrupole-quadrupole and the cross terms 
dipole-quadrupole. The dipole-dipole and quadrupole- 
quadrupole terms will be the same for atoms 1 and 2 
because they involve only the squares of the wave func- 
tions and therefore the sign difference in the p component 
for atoms 1 and 2 disappears. Therefore we concentrate 
only on the dipole-quadrupole cross terms which for each 
of the two atoms is given by: 
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FIG. 2: Ge partial density of states and the intensity of 600 
reflection. The DOS and structure factor are broadened by 3 
eV Lorentzian. 



The two terms are obviously not equal yielding a net 
result for each. For the K-edge we see that the r^ terms 
get us from the ground state to the d component of 
the excited state eigenfunction and the dipole (r) term 
gets us back from this eigenstate via the p component 
to the ground stat or visa versa. Immediately we see 
that the phase of the p-d mixing in the eigenstate en- 
ters and therefore the two Ge atoms discussed above will 
in fact scatter differently. Now one might say that this 
would also be the case for the s-p hybridization but a s-s 
transition is not quadrupole allowed because an s orbital 
has no quadrupole moment. This kind of hybridization 
could only be observed in a combination of a dipole and 
monopole transition. 

Now that we understand which physical scattering 
mechanism can distinguish between these two Ge atoms 
we can check if there is such a difference seen in the elec- 
tron structure of Ge. To do this we look at the conduc- 
tion band part of the band structure in a LDA calculation 
(TBLMTO-47 computer code (l3|), where we neglect the 
core hole potential. The 4p projected local density of 
states would correspond to the energy dependence of the 
K-edge x-ray absorption spectrum broadened by the core 
hole life time. This core hole life time is about 3 eV. The 
600 and 222 reflection structure factors in this mecha- 
nism will be determined by the energy dependence of the 
product of the amplitude of the p and d component of 
the band eigenstate times a radial integral and is given 
by: 
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where G = K' — K, m and m' are magnetic quantum 




FIG. 3: The geometry of x-ray scattering, 
number of d and p orbital respectively 
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j is a band index, C-I^{k) are components of the eigenvec- 
tor projected on a particular orbital of site n with lattice 
vector Rn and Yim are cubic spherical harmonics. 

In Figjl we show the Ap projected partial density of 
states and this convoluted with a Lorentzian of width 
3eV. We note that we have not taken into account higher 
lying p bands than the 4p band to which dipole tran- 
sitions are of course also allowed so the experimental 
absorption edge will fall off somewhat less rapidly than 
calculated here. To calculate the resonant amplitude we 
use the scattering geometry of Fig.^ and the angles deter- 
mined by the Bragg conditions for the given energies. At 
a photon energy equal to 11107 eV the scattering angle is 
36.35 degrees for the 600 reflection and 20.01 degrees for 
the 222 reflection. The resonant intensity will be deter- 
mined by the product of the p and d amplitude at a given 
energy in the band structure integrated over k space and 
summed over atoms 1 and 2 with the appropriate phase 
factor e"**^'^". As discussed above the product of the 
p and d amplitude for atoms 1 and 2 will have opposite 
sign which in the sum will be cancelled by the opposite 
phase factor for the reciprocal wave vector 600 or 222. 

In Figure g we also show the 4d partial density of states 
and we see that in general the Ad states are as expected 
higher in energy than the Ap states. It is then easy to see 
that there will be a phase change in the p-d mixing as we 
go from the edge to the high energy states at about 30 
eV higher. The doted dashed curve is then the resonant 
scattering intensity calculated for the scattering geome- 
try of Fig.H which demonstrates the strong energy depen- 
dence of the p-d mixing resulting in a strong peak close 
to threshold and a weaker structure 16 to 18 eV higher 
in energy and another peak at 30 eV above threshold. In 
fact the energy scans done to date did not cover the high 
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FIG. 4: The energy dependence of the phase of the 600 reflec- 
tion in a-a (tt-tt) channel with < < 90 or 180 < < 270; 
CT-TT with -40 < (^ < 40 or 140 < <;/> < 220; tt-ct with 
40 < <^ < 140 or 220 < (j) < 320. To obtain the phase for 
other (j) angles one needs to shift these by 180 degrees. 



energy region so this should be considered as a prediction 
which should be checked experimentally. The first strong 
peak at threshold agrees very well in intensity and width 
with the experimental data of Lee ei a/ |ll| . 

We notice that the resonant intensity goes to zero at 
about 1135eV which is a strong indication that the phase 
of the mixing changes sign. In order to study this we show 
in Fig.y which shows a strong energy dependence and a 
change of sign at about 18 eV above the edge and again 
at about 30 eV above the edge. In fact we predict another 
resonance at higher energies and this should be looked for 
in the future. In addition to the prediction of the second 
resonance we also predict that the phase of the scattering 
changes by 37r. This change of phase should change the 
interference with the omweg scattering [n4| dramatically 
and also this should be studied as a function of energy. 

In agreement with the paper by Templeton and Tem- 
pleton PI our calculations show that only the dipole- 
quadrupole term contributes to the resonant intensity in 
the crystal structure of diamond. First of all, the pure 
dipole as well as quadrupole term in Eq.(Q|) is essentially 
proportional to local partial p or d density of states pro- 
jected on a particular site. Due to the fact that there is no 
difference in the correspondent partial densities of states 
of the 8 germanium atoms in the unit cell and the reverse 
of signs in structure factor the intensity of resonant scat- 
tering in dipole-dipole or quadrupole-quadrupole chan- 
nels is zero. Secondly, as already mentioned, the mixed 
dipole-quadrupole term involving transition from Is to 
Ap and from 4s to Is orbitals of germanium also docs not 
give any contribution to the resonant scattering because 
neither of these are quadrupole allowed. This means that 
sp3 hybridization will not contribute to this scattering. 

The net tensor has then only two significant contri- 
butions which, in for example the p-d dipole-quadrupole 
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FIG. 5: The azimuthal angle dependence of resonance peak 
at energy 11109.2 eV. The upper panel shows 600 reflection 
for scattering angle 9 — 36.35°. The lower panel shows 222 
reflection at 61 = 20.01°. 



transition has a symmetry: 

'0 c o^ 

C 

^0 C Oy 

Here C is some constant for each particular energy 
and independent on azimuthal angle (f). The rows and 
columns of this matrix are numerated according to quan- 
tum number order (m — —1...1). The reverse transition 
has exactly the same symmetry and sign but both of them 
have opposite sign to the matrix elements related to the 
second type of germanium atoms in the unit cell. The 
latter is a direct consequence of wave function symmetry 
caused by the particular mixing of p and d orbitals in the 
diamond crystal structure. 

Finally we present the calculated azimuthal angle de- 
pendence of the 600 and 222 reflections at K-edge. Fig.^ 
shows that the intensity of the 600 reflection follows a 
(sinip cosip)'^ behaviour in a-a and tt-tt channels with 
equal to zero when the (001) axis is in the scattering plane 
and the incident beam is approximately parallel to this 
axis Fig.||.For the 222 reflection it changes to more like 
cos^cj) behaviour whereas the period in ct-tt (tt-ct) channel 
is changing from tt to 27r. 

For the 600 reflection we also calculate the absolute 
value of the structure factor following the paper by 
Colella and Merlini ll^ where we choose the polariza- 
tion factor and the Debye- Waller factors to be 1. To 
compare our calculations with experiment we define the 



integrated intensity / as ^^ I^-n + ^a-a were s is a spin 
index. In that case for an azimuthal angle of 25 degrees 
the integrated intensity is calculated to be 1.656 • 10~^° 
and the structure factor 0.44 which is in excellent agree- 
ment with 0.475 determined from the experimental data 
by Colella §. 

In conclusion we have demonstrated that the 
quadrupolar-dipolar resonant scattering in Ge will in- 
deed yield a 600 scattering intensity and that this is a 
direct measure of the p-d hybridization as a function of 
energy in the conduction band. Band theory docs an 
excellent job of explaining both the magnitude and the 
energy dependence of the resonance. In addition we have 
made a number of predictions regarding more resonances 
at somewhat higher energies and we have predicted the 
phase. Future experiments related to these predictions 
should provide an even more critical test of the theory. 
Since such p-d mixing is only allowed if the atom is not 
at an inversion center this mechanism could also be im- 
portant for all such materials and may explain part of 
the structure seen at resonance for the transition metal 
oxides in general and V2O3 in particular. 
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